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. We investigate a mean-field approach to a quantum brownian particle interacting with a quan- 

^ ' turn thermal bath at temperature T, and subjected to a non- linear potential. An exact, partially 

, classical description of quantum brownian motion is proposed, which uses negative probabilities in 

■ its intermediate steps. It is shown that properties of the quantum particle can be mapped to those 
of two classical brownian particles in a common potential, where one of them interacts with the 
quantum bath, whereas another one interacts with a classical bath at zero temperature. Due to 
damping the system allows a unique and non-singular classical limit at ?i — » 0. For high T the sta- 
tionary state becomes explicitly classical. The low-temperature case is studied through an effective 

■ Fokker-Planck equation. Non-trivial purely quantum correlation effects between the two particles 
are found. 

S 

> 

^ ' The main conceptual problem in quantum mechanics remains the link between the quantum and the classical 
I , worlds. Therefore, significant efforts were made over years to understand at least part of the quantum world in 
' ^ ' classical notions. Curiously, the quasiclassical domain, which should be the main subject of this understanding, is still 
itself under extensive investigation. Indeed, it is known to be non-trivial; in a sense it can be even more complex that 
^ ' the classical and quantum extremes alone. It is important to realize in this context that the quasiclassical domain is 
, not exhausted by the conventional ansatz h —> p]-p|, since this limit is singular (therefore some a priori concepts 
similar to coarse-graining are sometimes involved B ) , and since it does not commute with other limits of physical 
^ interest, e.g. the limit of large times. 

One of the established approaches to the quasiclassical domain is a collection of mean-field methods known as 
[ Gaussian decoupling procedure or time-dependent variational approximation Q-j^. Mean-field (variational, Hartree- 
• Fock) methods are well-known in the quantum theory, and were applied for a while in many different areas. This set 
of methods appeared to be especially suitable for the quasiclassical domain, since it attempts to realize in a simple 
and straightforward way the above-mentioned program of understanding the quantum theory in classical terms. An 
' impressive amount of experimental confirmations, in particular in quantum chemistry and atomic physics Q , numerical 

■ and self-consistency checks were made for those methods. Therefore, they have already become a well-formulated and 
sound approach. 

^ ■ In the present paper we apply this mean-field method to the simplest quantum dissipative system: A quantum 
r ' brownian particle interacting with a thermal bath. There are specific reasons to study this type of quantum systems in 
^ the context of the above-mentioned problems. In contrast to closed Hamiltonian systems, quantum dissipative systems 
Q [ are more reliable candidates to understand their physics in classical terms. Indeed, a non-unitary evolution provides 
O • a natural mechanism of decoherence [l^ , and the need of artificial coarse-graining procedures is gotten rid of. In the 
J> I light of this conceptual advantage it should be surprising that the basic understanding of their theory is still rather 
fragmentary. This is so mainly because one can use neither general properties of unitary evolution, which describes 
. , closed systems, nor markovian properties of the classical stochastic dynamics, since due to the correlation time h/T, 
I which is relevant at low temperatures, the corresponding statistical dynamics is essentially non-markovian. Both 
these dynamical properties provide important general information, which is, thus, absent in the quantum dissipative 
case 1^. The only exceptions are weakly damped high-temperature systems, where a copulation of those two things is 
possible, namely, an influence of the bath is described classically, whereas the rest remains quantum-mechanical [ p^ . 
We will not be concerned with them within the present paper. A fairly general theory of strongly-damped and/or 
low-temperature dissipative systems is still under construction, though some suggestive results were obtained recently 
[ p^ , where for a class of systems two of us proposed a consistent statistical thermodynamical theory of quantum 
brownian motion. 
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Our plan is the following. In section II we will first briefly recall some known facts on quantum Langevin equations 
and quantum noise. Here we obtain also generalized Wigner-Moyal and von Neumann equations, which are exact 
consequences of the quantum dynamics, but describe the brownian particle in almost classical terms. The Gaussian 
approximation will be presented in section III, where we will obtain the basic equations of the present paper, and 
draw some general conclusions. Here we will discuss, in particular, what is the constructive role of friction when 
establishing quantum-classical transition. The effective Fokker-Planck dynamics will be discussed in section IV. In 
section V we study low-temperature properties of the model. We conclude in the last section. 



II. QUANTUM LANGEVIN EQUATION AND LIOUVILLE-MOYAL EQUATION 

A. Quantum Langevin equation 

This fundamental equation of quantum brownian motion theory is derived from the exact hamiltonian description 
of a subsystem (brownian particle) and a thermal bath, when tracing out the degrees of freedom of the bath. The 
standard assumption is that at the moment t — the states of the subsystem and the bath were decoupled from each 
other, and the bath was in equilibrium at temperature T |p|,p|Jl4|. Further, the influence of the particle to the bath is 
assumed to be sufficiently small; thus, only the linear modes of the bath are excited, and the interaction between the 
particle and the bath is linear. Since the dynamics of the bath is linear, it can be solved exactly. Following this line 
of exact calculations , one derives the quantum Langevin equation 

P 

a; = — , 
m 
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m 



p + — dt'-t{t~t')p{t') + v'{x)^-jre-'^'x{o) + 'n{t), (1) 







where p{t) x{t) are Heisenberg operators of momentum and coordinate, and V{x) is an external potential. The 
parameter 7 is the damping constant, which determines the interaction between the bath and the particle. For 7 — s- 
one gets from Eq. (Q) the usual Heisenberg equations. F is the maximal characteristic frequency of the bath, and it 
determines the retardation time of the friction kernel 

7W=7re-r|*l- (2) 

The operator fjit) is the random noise, which appeared due to the uncertain character of the initial (equilibrium) 
distribution of the bath. This noise can be shown to be gaussian, due to the fact that the thermal bath is a harmonic 
system and was in equilibrium. It has the following properties: 

Kit) = l{mm+mm)n ^ l{m;m)^ = ^£ ducocoth cos{u:t)^-^^-^, (s) 

mm - mm ^ m),m] - (4) 



Hereafter we ase 



{A;B) = ]^{AI3 + I3A), [A,B] = AI3 - BA, [A,I3]+ = AI3 + BA (5) 



for any operators A, 13. 

The connection between properties of the noise and the friction kernel is the consequence of quantum fluctuation- 
dissipation theorem ||^ . Eq. (|l|) with physically suitable forms of the potential and friction describes a rich variety of 
physical phenomena (see references in p|,p|,p^ ) . 



1. Quasi- Ohmic limit 

In the present paper we shall restrict ourselves to the quasi-Ohmic case, where F is much larger than other char- 
acteristic times, but still finite. The main reason of this approximation is to have an exact equation for the Wigner 
function of the brownian particle, which will be derived in the next subsection. 
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For the quantum noise one has that K{t) = — In(rt) > for small times, and for t ^ l/F, K{t) is anticorrelated 
with the universal correlation time h(3/2TT: 



K{t) = 



sinh 



(6) 



Being coherent, the low-temperature quantum thermal bath neccesarily generates a colored noise. The classical white 
noise situation is recovered when taking the high-temperature limit (hf3 ^ 0) . In general, this should be done before 
the limit T oo. Notice that in contrast with the classical case, the quantum noise does not disappear for T ^ 0, 
since even in this limit the initial state of the quantum thermal bath remains indeterminate. 
In the quasi-Ohmic regime one can expand the memory kernel for the friction in Eq. ([^) 

di'e-^(*-*')p(i') = - fe-^'PiO) + ^e-^'m. (7) 

For t > and large F the exponentially small factors depending on the initial conditions can be omitted in Eqs. (|^, 
0), and we finally get 

p+^p{t) + V'{x)^^it), (8) 
m 

In other words, for not very small times one can regularize only the noise, but keep the friction local, when considering 
the quasi Ohmic case. 

Let us finally notice the following general comutation relation between the noise and an arbitrary operator A(t) of 
the particle 0,1: 

[Ait),fj{s)]^ j\t'[Ait),xit')]^^lis-t') (9) 

Here the coordinate operator appeared just because the interaction with the thermal bath is taking place through the 
coordinate 0,1). In the quasi-Ohmic case this will go to 

[A{t), fiis)] = 27^{^(i - s)[A{t), x{s)]}, (10) 

where 6{t — s) is defined to be 1/2 at t = s. It is seen that [fi{t),x{t)] = [f)(i),p(<)] = 0. One could attribute these 
relations to causality, but they are an emergent property of the quasi-Ohmic limit. 



B. Generalized Wigner-Moyal equation 

The quantum Langevin equation (|^) is a non-linear operator equation, and as such it can hardly be handled directly. 
Here we will present a generalized Wigner-Moyal equation for the corresponding Wigner function of the particle, which 
exactly corresponds to Eq. (^) in the same sense as the pure Heisenberg equations correspond to the usual Wigner- 
Moyal equation. Besides technical advantages which will be used further, this equation presents an interesting account 
for an exact description of the non-linear quantum problem through proper classical terms. 

One is looking for an equation for the Wigner function: 

W{x,p,t) = (trpo IF(x,p,t))^, 
^ f da db 

Wix,p, t)-j^ eM-^ax - ^bp + ^ax(t) + ^bpm (H) 

where po is the initial (at the moment t = 0) density matrix of the brownian particle. W{x,p, t) is the support of the 
Wigner function, which can be viewed as a quantum analogue of the delta- function. In the classical limit, where x and 
p approxmately commute, W{x,p^t) tends to the ordinary probability distribution of the coordinate and momentum. 

Since the equation for W{x^p) is expected to be linear, one is interested by an equation for W{x,p), whereas the 
averages can be taken later. The derivation of this equation is fairly straightforward, since it uses Eqs. (||, 11) and 
the standard commutation relation [x{t),p{t)\ = iTi. We will write only the final result: 
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n— 1 ^ ^ 

The first three terms in the r.h.s of tliis equation are the standard drift terms of the classical Liouville equation, and 
the sum represents a purely quantum correction, which comes from the non-linearity of the potential. 
Our object of interest is the last term, which upon averaging will look like 

d_ 

dp ' 



— {^{t)W{t))f,. (13) 



Notice that ?7(i), W{t) commute, since fi{t) commutes with x(t), p{t), and W(t) can be presented as a product of two 
terms, which depend only on x{t) and p{t) correspondingly. 
Let us adopt the following formal expansion for W{t) 

W{t\Tj)^W{t\0) + -y2 dsi...ds„i?(si,..,s„)n[?7(si)...rKs„)], (14) 

n[??(si)...??(s„)] = — 77(s,J...77(si„). (15) 

Here n[77(si)...77(s„)] is the symmetrized product; the coefficients i?(si,..,s„) are c-numbers. Due to the fact that 
[W{t),ri{t)] = one has: 

{fi{t)W(t))f, = {fi{t);W(t))f, = Y,- / dsi...ds„i?(si,..,s„)(^(i);n[7)(si)...7)(s„)])^. (16) 



n! 

n=l 



Since f)(i) is a gaussian random operator with {f]{t))fi = 0, one can use Wick's theorem: The correlation of an odd 
number of ffs vanishes. The correlation of an even number of ffs is equal to the sum of products of pair correlations, 
the sum being taken over all pairings. For example: 

{f]{ti)f^{t2)ii{t:i)7^{u))fi = mi)f]{t2))fi{f]{h)Ti{u))f, + mi)vih))f,m2)f,iu))^ + mi)fi{u))f,m2)f,{h))v (17) 

In this way one derives: 

n 

(??(i);n[?7(si)...r)(s„)])^ = Y{fl{t);fi{sa))fj {n[fl{si)...fi{sa-i)fi{sa+i)---fi{sn)])fi- (18) 



Having substituted this equation to Eq. (16) one gets 



{fi{t)W{t))fj ^y^l TTT / dsi(f)(t);f)(si))^ / ds2...ds„i?(si, .., s„)(n[7y(s2)-?7(sn)])^. 

„^il"-tj-Jo Jo 



(19) 



Now one notices that the only feature of the quantum noise which enters here is the autocorrelation function 
{'rj{t);'fj{s))fi = Kit — s) because in the end (Jl[r){s2)--.'q{sn)])f, can be expressed through it using Eq. (^ several 
times. So nothing will change if we replace f}{t) in Eq. ( p^ ) by a classical gaussian noise rj{t) which has the same 
autocorrelation function K{t — s), 

mw{t))f, = mw{t))^. (20) 

Now we substitute this result into Eqs. ([l^), take the average over the initial state, but take out the average over the 
classical noise: 

71 — 1 ^ ' 

where the true Wigner function I^(x,p, t) will be obtained by averaging over the classical noise 

W{x,p,t)^{w{x,p,t))^, (22) 

where w{x,p,t) is an auxiliary object. Thus in the quasi-Ohmic limit one can treat the quantum noise as a purely 
classical object as far as the Wigner function is concerned. Notice that this result is exact and relies only on the 
quasi-Ohmic limit. In that respect it is different from the description through semiclassical Langevin equations, where 
in the overdamped (large 7) limit one also gets an analogue of Eq. (^) but without the ?i-dependent terms. 
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C. Generalized von Neumann equation 



Here we will investigate the von Neumann equation for the density matrix, which corresponds to the Wigner function 
w{x,p): 



Since Eq. ( |2l| ) is linear, one directly obtains for r{x,x',t) = {x'\f{t)\x) 



ih—r(x,x ,t) 



rix,x\t), (24) 



or in the equivalent operator notations: 



— ^--[H~x7^,r] + ^[{p,f},x]. (25) 

To obtain this equation from Eq. ( pT] ) in quick way, one can use the following correspondence between operators A, 
B and their phase-space representations A^'^\x,p)^ i?^^'(x,p): 

= dpA^"^ d^B'^"^ - dpB^"'^ d:,A^''\ [A,B]+ = 2A(")b('=), (26) 

and take into account that x^'^^ — x, p^'^-' = p, r*^^^ = w{x^p). 

The true density matrix /) can be obtained after averaging by the classical noise 'q{t): 

m = m),- (27) 

Therefore, all possible averages are obtained as 

tr(pi) = (tr(ri))^, (28) 

where an operator A lives in the Hilbert space of the brownian particle. 

It should be stressed that f is not a density matrix itself, but rather a tool to calculate averages. Indeed starting 



from Eq. (25) one easily gets 

— trf = ^ trf(i) = 1, (29) 
at 

Atr(r(t)2) = -trifitf) ^ tr{f{t)^) = e^*/" tr(f(0)2). (30) 
di TO 

These indicate that r has negative eigenvalues for i > 0, the absolute value of which grows with time. However, 



positive eigenvalues compensate this growth in such a way that Eq. (|29| ) holds. In particular, Eqs. (29, 30|) show that 
a dissipative systems cannot be described by a wave function even if the averaging over the noise is postponed. 

Thus we see that the "classicalization" of the noise is not just a technical procedure, but it has to be accompanied 
with a change of interpretation: The unaveraged density matrix f is not a true density matrix, since it does have 
negative eigenvalues. In other words, explicitly classical components of the dynamics lead to the appearance of 
negative probabilities. On the other hand, there are no reasons to consider those negative eigenvalues as something 



unphysical: Our derivation of Eqs. (^ 25) was exact, and later we will present other indications that neither the 
true density matrix, nor averages calculated according to Eq. ( |2^ ) show unphysical properties. The situation is the 
same as for the Wigner function at a given time, which is not a positive probability density, but can be safely used to 
evaluate expectation values by integration. 

The appearance of negative probabilities as a result of imposing partially classical properties has certain analogies 
with the classical interoperation of quantum entanglement, and in particular Einstein-Podolsky-Rosen phenomenon. 
This interpretation also uses negative (though not directly observable) probabilities [TT| . 

Finally we would like to mention that the reported non-positive character of the unaveraged density matrix has 
nothing to do with the known technical problem which also appears through non-positive (averaged) density matrices 
in certain quantum markovian diffusion equations [ p2[ . There the problem is merely technical and arises due to the 
fact that the markovian approximation in the quantum theory of open systems is essentially time-inhomogeneous, 
so that its careless use leads to such problems. Indeed, the problem disappears after a more consistent treatment of 
the situation p^ . In our case, in contrast, the non-positive character of the unaveraged density matrix is an exact 
consequence of our attempt to handle the quantum noise classically. 
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III. GAUSSIAN DECOUPLING PROCEDURE 



Eq. ( pTl ) is, of course, intractable in general. It joins all technical difRculties of the classical Liuoville equation with 
a given noise and friction and those of the pure Moyal equation. In other words some substantial simplifications are 
necessary to proceed further. Here we will apply the mean-field approach, namely, a solution of Eq. (21 ) will be looked 
through its moments, and the gaussian decoupling procedure will be applied to the higher-order moments. By its 
spirit this is very similar to the Grad method in the kinetic theory of rarefied gases [|l6| , and has been applied recently 
in quantum theory of closed systems as well [^|j4j. It is clear that the consistency of this approximation should be 
checked together with the final solution. At the moment we will notice only that since its application is connected 
with the weakness of quantum fluctuations, it will have a reliable range of validity in the quasiclassical domain. 

The working variables will be 

(31) 
(32) 
(33) 
(34) 
(35) 



dx ~ (x) = J dp dx xw{p, x), 

dp = {p) = J dp dx pw{p, x), 

dxx = {{x - {x) f) = j dp dx [x- {x))'^w{p,x), 

dpp = {{p - {p)f) = J dpdx {x~ {x)fw{p,x), 
1 

2^ 



dxp = T^ii^ - {x)){p- {p)) + {p- {p)){x - {x))) = I dpdx {x- {x)){p- (p))w{p,x). 



For the higher-order correlations one assumes the gaussian decoupling: 



dp dx w{p, x){x ~ {x)f" = {2n - 1)11 



dp dx w{p,x){x — {x)Y 



dp dx w{p, x){p — (p)){x — (x)) 



2n+l 



(2n + 1)!! 



dp da; w{p, x){p — {p)) 



dp dx w{p,x){x — {x)y 



(36) 



(37) 



This just means that w{p,x) is restricted to the subspace of gaussian functions: 

1 



w {x,p\dxx,dxp,dpp,dx,dp) = ^ exp [dpp{x ~ dx)^ + dxx{p - d^f - 2dxp{x - dx){p - dp)]^ 



dxx^pp d'xp' 



(38) 

Applying Eqs. (|31|-|37D in Eq. ( ^l| ) one gets the following equations of motions, which are now classical equations for 
classical variables: 

d , dp 
dt m 



_dp = -^dp+vit)-v\dx) + j2 



1 

m 



V^'-+'Hdx) 



(2n)\\ 



at m 

dt ~ m ^ (2n - 2)!! 

n— 1 ^ ^ 



d _7 1 

rif m m < ^ 



v^'^Hdx) 



dr-^'^ m""'' ' m"'^'^ ^(2n-2)!! 

n— 1 ^ ^ 

This set of equations can be considerably simplified if the following change of variables will be made: 

dx — X, 
d„ = P, 



(39) 



(40) 
(41) 
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t2 2 

n a 



(42) 
(43) 

dxp = Qn, (44) 



dpp = ^+^^ (43) 



where a is chosen such that 





P 

TO ' 






-^P + 7y(t) - 
m 


" ax' 




n 

m' 






TO i9Q ' 




d 


TO 


= e-''*/"cr(0) 



i = tr(r2) (45) 

(7 

With this change of variables Eqs. (^9|) will read 

(46) 
(47) 
(48) 
(49) 

i-"^*/", (50) 
where the initial state was chosen to be pure for simplicity: (7(0) ~ 1, and where 7i is an effective Hamiltonian: 

n— 1 • ^ 

The true Wigner function for the original quantum particle will read: 

W{x,p)^ J dX dQ dP dU w{x,p\X,Q,P,U)V{X,Q,P,U), (52) 

where 'P{X, Q, P, U) is the ordinary probability distribution of the classical random variables X, Q, P and 11. 

The physical meaning of this approach is now clear. The Hamiltonian ( pl| ) corresponds to two classical particles 
with coordinates X, Q and momenta P, H. Eqs. (|4^p9|) show that X-particle couples to the quantum bath through the 
damping jP/m and noise r]{t). Although r](t) is not an operator, its correlator is still given by the quantum spectrum 
K{t). Q-particle interacts with a classical bath at zero temperature, since only in this case a classical particle is 
subjected to damping, but not to noise. The effective Hamiltonian Ti. non-trivially couples these two particles. It is 
also time-dependent, though this dependence is quite simple. 

Already the general form of the effective Hamiltonian leads us to the following important observation. As follows 
from the derivation of the result, the purely classical case corresponds to ?i ^ 0, H ^ 0, q — > 0. Without damping 
one will have a — const, which just reflects unitary evolution, where is conserved. Since Q, H typically have order 
0{h), all terms in Eq. (^) should disappear in the classical limit. For all terms besides the last one this disappearance 
is clear. This last term can make the classical limit non-unique or even singular. This phenomenon is well-known in 
quantum chaos Moreover, even fairly simple integrable systems can display singularities in the classical limit 

Now we observe that this dangerous term h'^a{t)'^ / {8mQ^) disappears with the characteristic time 7/to, thereby 
ensuring the relatively straightforward classical limit in a damped system. Notice in this context that usually it is 



only the noise, which is believed to facilitate the classical limit, providing a mechanism for decoherence |13|. 
Notice that the equality 

dxxdpp — d^p = - , (53) 

when a given by Eq. ( pO[ ) is less than one, does not indicate a breaking of the uncertainty relations, since f itself is 
not a density matrix. The correct uncertainty relation will read as 
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where ((...)),, is the complete average, namely the average by r (indicated with (...)), and by classical noise (indicated 
with (...),)). 

These observable averages can be expressed as 

((x2))„ - {{x))l = (d,,), + {dDr, - {d.,)l (55) 

{{P^))v - ml = {dppU + {4)v - {dp)l (56) 

^((((^ - {{x))v){P - {{P))v) + iP- {{P))v)i^ - {{i))ri)))v) - {dxp)r, + {d^dp)^ ~ {d^)r,{dp)n (57) 

Recall the situation with the exactly solvable harmonic potential. Here dxx, d^p, dpp tend to zero in the long-time 
limit, being decoupled from dx, dp. Then Eq. (^) is obviously satisfied. 

Let us briefly mention another aspect of the proposed scheme, which can be interesting on general grounds. Two 
important length-scales are associated with any quantum system 



Lc = VD^, (58) 

= "l^iDpfD^I-DlpY ^^^^ 

where D^x, Dpp are dispersions of the coordinate and momentum, and Dxp is the corresponding cross-correlation. 
These lengths quantify the quantum {Lq) and classical (Lc) aspects of the system, since due to uncertainty relation: 
Lc/Lq > 1, and the classical limit corresponds to Lc/Lq 3> 1. 
The analogous lengths for our variational scheme read 



n'dx:. 



LjLq ~ e-^*/". (61) 

The lengths ic, Lq still characterize classical and quantum effects, but the uncertainty relation does not apply to 
them, due the above remarks. Ratio (^) can now be much smaller than one, namely the quantum effects can be 
overdominating. The noise is needed to recover the uncertainty relation and to limit the overspread of quantum 
effects. 

To conclude this section we will notice that recently an attempt was made to construct a mean-field theory for a 
quantum brownian particle [p^ . However, the authors did not start from the correct statement of the problem, and 
were led to an incorrect result that their open quantum system can be still described by a wave function (i.e., a pure 
state). In our notations this will amount to put a{t) = 1 for alH > which is clearly incorrect in the light of Eqs. ( |5C| , 
|30| ). Besides the technical aspect it contradicts to the general qualitative statement which we draw after Eq. (|50|). 

IV. FOKKER-PLANCK EQUATION 

Eqs. ( ^6[]49| ) are still fairly complicated non-linear equations. To study them especially at low temperatures we 
shall employ methods recently developed by two of us jl^ . 

We are looking for an equation describing the common probability distribution 

^(2/1, 2/2, 2/3, 2/4, i) = {Sivi ~ P{t))5{V2 - X{t))5{v:i - n{t))5{vi ~ Q{t))),, (62) 
Using Eqs. ( ^6[]49| ) and direct differentiation one will get 

dV d{vkP) d^s{P{t) - y,)S{X{t) y2)S{Uit) - ys)5{Q{t) - y,)v{t))„ (63) 



dt ^ dyk dyi 



where 
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m 

yi 

V2 = , 

TO 

ys 

V4 = . 

TO 



Since the gaussian noise is distributed with a functional 

1 



n[7^] - exp-| / dtdsr^{t)K-\t-s)7j{s), 



one has 



dt'K{t' - 1) 



'Sriit'Y 

Substituting this equation into Eq. (|6^) one obtains after functional integration by parts: 

4 

dt " 



^ dyk dyi \ ST]{t') 



- {diP{t) - yi)S{X{t) - y2)Sim - y3)S{Q{t) - ^4)} 



(64) 
(65) 
(66) 
(67) 

(68) 

(69) 

(70) 



To calculate the functional derivatives entering this equation, we notice with the direct variation of the equations of 
motion: 



6P{t)/5v{t') 



6{t - t') I exp / duA{u) 



(71) 



A{x) = 



where {...}+ means the chronological ordering. 

Substituting the last expression to Eq. (R3) one gets 





— 7/m 







-dlQ-H \ 




1/to 


















—7/m 




V 








1/m 


0/ 



yi)5{X{t) - y2)5m) - y3)HQit) - yi)^ki{xm , 



where $ is the following 4x4 matrix 



dt'K{t') jexp^ d 



uA{x{u)) 



(72) 



(73) 



(74) 



and <i>fei is the corresponding matrix element. This result is still exact, but untractable, since it involves the functional 
$ of the history {X{t')}, {Q{t')} for t' < t. In the classical limit one gets for t' > the white noise K{t') 2^T5{t'), 

^ki = iTSki, 

thus reproducing the corresponding classical Fokker-Planck-Kramers equation. A closed equation for V can be ob- 
tained also in the harmonic case, where A does not depend on X(u), Q{u). These two exact realizations prompt the 
way to proceed in the nonlinear case. Since Kit') exponentially decreases for t > max(ft,/3, r~^), and this time can 
assumed to be small in the quasiclassical domain, one can make a Taylor expansion of the exponent in Eq. (|7^), and 
keep only the first term: 
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exp 



duA{x{u)) 



(75) 



Due to the 5-function in eq. (|73|), we may then replace the fluctuating x{t), q{t) by the sure variables a;, q, after which 
$ is no longer a fluctuating quantity, and can be taken outside the averaging in eq. (^3|), thus bringing a closed 
equation for P. $ will be calculated with help of the folowing formula 



sA{x) 



kl 



dz 
2m 



1 



A - z X 1 



(76) 



kl 



where 1 is the 4x4 unit matrix. 

Our step ( |75| ) is still exact for the harmonic initial potential V , while in the general case with a characteristic scale 
of anharmonicity L a condition 



»1, 



(77) 



is to be satisfied, which restricts the correlation time hf3. If this time is not small enough, one can notice that the 
linear part of A already suppresses exponentially the large values of s, and make the same expansion in Eq. (^3|). This 
is valid when the nonlinearity of the potential is small with respect to its linear part. Notice a certain nonperturbative 
aspect of the result, since the expansion was made inside of the exponent. 
The final result that we obtain is a diffusion-type equation for P itself: 



dV{P,X,Tl,Q,t) 



dt 



-fDppiX,Q,t) 



P dV 
m dX 



d_ 

dP 



Ip- 



m 

dX 



m oq 



ap2 dpQx 



[Dxp{X, Q, t)V] + iDnp{X, Q, t) 



d_ 

dPdIi ^ dPdQ 



6)2 



m 

dQ 



V 



[DQp{X,Q,t)V], 



(78) 



where we have changed yi ^ P, y2 ^ X^ j/3 i-^ 11, j/4 i-^ Q, and diffusion coefficients Dpp, Dxp, Dup, Dqp are 
instantaneous functions oi X , Q and t, and no longer functionals of the history. Since the analytic structure of the 
diffusion coefficients is somewhat involved, it will be explained gradually. All diffusion coefficients converge to finite 
values for large times. This convergence is exponential, and has a characteristic frequency 



mini<fe<4[Re(a;fe)], 



(79) 



where 



1 ± 



4m&i 



(80) 



2 m 



i± Wi 



(81) 



and where 



(82) 



It seen that in order to have YieLOk > 0, which is necessary for convergence, one has to require the conditions of local 
stability 



dxxn + OqqU > 0, 
dxxH dxQ-H > [dxQ-Hf. 

Hereafter they will be assumed to be satisfied. 

Let us now present explicit formulas for the stationary values of the diffusion coefficients: 



(83) 
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Dpp{X,Q) 



1 



duj 



m{bi ~ 62) Jo Ti" 



bl + dggH 



dqqli. 



(84) 



Dxp{X,Q) 



1 



to(6i - 62) Jo 7^ 



(61 + a,,7i)(cj2 + 61/m) (62 + dqqH){Lo'^ + b2/m) 



(85) 



^np(^,Q) 



m{bi-b2)Jo TT [(Cj2 +tj2)(^2 _^^2) (tj2 + ^) (^2 ^ ^2) 



(86) 



m(6i - 62) Jo 
where K{uj) is the spectrum of K{t), 



up' + bxjm, 



a>2 -|- b^jm 



(c^2+^2)(^2+^2) (0,2 (^2 +^2) 



i^(.)=.?.7-coth(^— j 



(87) 



Let us indicate that the diffusion process ( [Tq ) is non-markovian. In the purely classical limit the time-dependence 
in the diffusion coefficients disappears, and all diffusion coefficients besides -Dpp T disappear as well. Only then 
Eq. ([tsI) describes a markovian process. 



V. REDUCED DESCRIPTION 

Since Eq. ( [ZS] ) is still rather complicated, it is reasonable to look for relatively simple limits. One of them is the 
overdamped limit, which is characterized by large values of 7. In this case Eq. ( [78| ) can be reduced to an equation, 
which describes only slow variables Q. To proceed with this limit we shall define the following moments 

MM{X,Q,t)= J dP dU P''U^ViX,P,Q,U,t), (89) 

and construct an equation for them starting from Eq. (|78|). 

1 1 (k + 1)^ 

Mki = dxMk+u - -dQMk.i+i - -Mk,i - kMk-u dxH - lMk,i-i BqH 

m m m 

+jk{k-l)DppMk-2j - kdx[DxpMk-u] - ldQ[DQpMkj-i]+lklDnpMk-u-i. (90) 
Let us write down few first members of this hierarchy: 

Moo^- — dxMio-—dQMoi, (91) 
m m 

Mio = --dxM2o - -dgMn - -A^io - MoodxH - 9x[^xpMoo], (92) 
m m m 

1 1 7 

Moi = dxMn aQA^02 Moi - MoodgH - dglDQpMoo], (93) 

m m m 

M20 = -—dxMao - —dQM2i - —M20 - ^MwdxU + 2-/DppMoo - 2dx[DxpMw], (94) 
m m m 

M02 = --dxMi2 - -9qMo3 - —M02 - ^MoidgH - 2dQ[DQpMoi] (95) 
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Mil = dxM2i OqMu Mil - MoidxH - MiodgH + -fDnpMao - dx[DxpMoi] ~ dqlDQpMw]. 

m m m 

(96) 

In the first order of large 7 one can skip the time-derivatives in Eqs. (|9l|, |9|), since they have al least order 0(7^^). 
Further from Eqs. ( p3[^3| ) one gets the following approximate relations 

M20 = mDppMoo, (97) 



X02 =0(7-'), 

Mil = ^mDnpMoo- (99) 

Notice from Eq. (|8^) that Dup is of order I/7 for large 7. These equations are substituted in Eqs. ( ^ , p2| ), which in 
combination with Eq. (^Tj) brings the following reduced equation for the Mqq = J-{X,Q,t), which is the probability 
distribution of the slow variables: 

^dtHX,Q,t) = dx[Tdxn]+dQ[TdQn]+dxQ[TDnp]+dxx[:FDxx]+dQQ[TDQp], (100) 
Dxx = DxP + Dpp. (101) 



The first two terms in the r.h.s. of Eq. (100) are due to drift, whereas other terms are responsible for the diffusion. 
Let us now discuss this situation in details. In the classical limit, which is realized for sufficiently large temperatures 
or for ?i — > 0, one has Dxp — > 0, Dqp — > 0, Dup and Dpp T. Thus, Eq. ( |100D goes to the corresponding 



classical Fokker-Planck equation. As Eqs. (|49|, |50|, 100) show, no noise is acting on the Q-particle, therefore in the 
classical case it just relaxes to zero and does not fluctuate at all. Therefore, despite non-linearity of the potential V{x), 
the classical variables X,P decouples from the quantum variables Q,Il and tend to the classical Gibbs distribution. 
The quantum variables disappear, as seen also from Eqs. (p8| , p9| ). Recall that the very fact of this homogeneous 
disappearance is connected with the exponential damping (BOp of the singular term a{ty / {8mQ'^) in the effective 
Hamiltonian (|5l|). 

On the other hand, in the quantum case Eq. ( |100| ) shows that both variables X and Q become correlated and 
involved in a common dynamics. Moreover, as follows from Eq. (p9|), there is a well-defined correlation between 11 and 
P. In particular, through interaction with the classical variables quantum variables become coupled to the thermal 
bath. 



VI. CONCLUSION 



This paper was devoted to the mean-field (variational) theory of quantum brownian motion. Mean-field methods 
are widely applied in quantum theory [^-^ and have an established range of validity. Their general property is 
reduction of an initially quantum problem to an approximation involving only effective (mean-field) variables with 
classical (commuting) behavior fl^ . The original quantum character of the problem is then refiected through an 
effective Hamiltonian. These properties of the mean-field description for closed systems were established also for more 
general cases (e.g. the time-dependent Hartree-Fock approximation), where the generated effective classical dynamics 
does not have the canonical form, but instead can be embodied into a more general Poisson structure fl^ . 

Being motivated by the effectiveness of the mean-field approach, we considered here its application to the problem 
of quantum brownian motion, which is the main representative model of quantum open systems. 

Our first step was to substitute the original quantum problem by an auxiliary semi-quantum one, where only a 
part of degrees of freedom is quantum. In this exact step it is possible to replace the operator- valued quantum noise 
by an auxiliary classical gaussian noise, which has the same spectrum as the original quantum noise. In that way we 
obtain the generalized Wigner-Moyal equation (|2^, |2^) . In this step our description uses negative probabilities in the 
sense that the unaveraged density matrix (^5|) does have negative eigenvalues. However, no unphysical results appear 
in the level of observables quantities. The negative probabilities appear as the cost for having explicitly classical 
elements in a quantum dynamics. This situation is reminiscent of the ordinary Wigner function which also cannot 
be interpreted as a probability density, but which shares some of its properties and does predict correct quantum 
mechanical averages. 
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At the second step of our description we applied the Gaussian approximation to the unaveraged Wigner function. 
By this procedure the initial quantum stochastic problem became reduced to a problem of two classical particles 
with friction and classical noise. This noise is nevertheless not white, but is correlated with the same spectrum 
as the original quantum noise. Further investigation allowed us to uncover an important role played by friction in 
establishing classical aspects of the problem. It appeared that a singular ?i-dependent term in the effective Hamiltonian 
( pT| ) is diminished by friction with a characteristic time (|o|) inversely proportional to the damping coefficient. This 
ensures the existence of the unique and homogeneous classical limit for times larger than the above characteristic 
time. This fact is contrasting with the non-commutation of the classical {h 0) and long-time t —^ oo limit for 
(closed) Hamiltonian systems Q. At low temperatures, where the quantum effects are essential, the dynamics of the 
mean-field degrees of freedom is essentially different, because they become correlated with each other. This is shown 



in particular by Eqs. (M, IOC). In other words, the effective classical dynamics still contains h/T as a correlation time 



of the classical noise, and therefore displays essential different dynamical behavior for high and low temperatures [ [14[ . 

In the present paper we restricted ourselves to the general framework of the mean-field quantum Brownian motion. 
More specific applications, e.g. for open many-body systems, are planned to consider in future. Finally, it is hoped 
that the paper will open a road for applications of mean-field methods in quantum dissipative systems. 
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